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Abstract 

In this paper, we prove via counterexamples that adding an advection term of 
the form Shear flow (whose streamlines are parallel to the direction of propaga- 
tion) to a reaction-diffusion equation will be an enough heterogeneity to spoil the 
increasing behavior of the KPP speed of propagation with respect to diffusion. The 
non-monotonicity of the speed with respect to diffusion will occur even when the 
reaction term and the diffusion matrices are considered homogeneous (do not depend 
on space variables). For the sake of completeness, we announce our results in a set- 
ting which allows domains with periodic perforations that may or may not be equal 
to the whole space R . 

Keywords: Traveling fronts, reaction-diffusion, monotonicity with respect to diffu- 
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1 Introduction and motivation 



Pulsating traveling fronts are particular solutions of heterogeneous Reaction- Advection- 
Diffusion equations that describe propagation phenomena for combustion models, evolution 
of epidemics, population dynamics and many other phenomena. This paper is dedicated 
to answer the question whether the minimal speed of propagation (KPP nonlinearity) is 
monotone with respect to the diffusion coefficient or not. We answer this question in the 
case where an advection term is present. This work is also a continuation of the paper 
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[I] which dealt with the dependence and asymptotic behaviors of the minimal speed of 
propagation on the coefficients of the reaction-diffusion equation in heterogenous media. 

To explain the question mathematically we start by recalling the definition of the min- 
imal speed of propagation and the changes that happen while passing from a homogeneous 
to a heterogeneous setting. 

1.0.1 Homogeneous framework 

In the article [6] by Kolmogorov, Petrovsky and Piskunov, the setting was "homogeneous" 
in the following sense. The equation considered in [6] was 

u t = Au + f{u) for all (t,x)eRxR N , (1.1) 

where < u = u(t,x) < 1 was a function defined over R x Mr. The reaction term / was 
assumed to be of KPP type. That is, / : [0, 1] i— > M such that 

f(u) > for all u G (0, 1), /'(0) > 0, (1.2) 

together with the sub-linearity condition (which is usually called the "KPP condition") 
over the interval [0, 1] 

VuG [0,1], f(u)<f'(0)u. (1.3) 

One can see that, in the equation (11.11) . the diffusion and reaction terms do not depend 
on space and time variables and there is no advection term of the form q(x) ■ Vw. Fixing 
a unitary direction e G M N , traveling fronts in the direction of — e and with a speed c, in 
such a homogeneous setting, are solutions u(t,x) = (p(x ■ e + ct) = 0(s) which satisfy the 
limiting conditions </>(— oo) = and 0(+oo) = 1. Kolmogorov et al. [6] proved that when 
the reaction term is of KPP type (i.e. satisfying the conditions (II. 2p and (11.31) ). there exists 
a minimal speed c* = 2a/ /'(0) such that a traveling front propagating in the direction of 
— e with a speed c exists if and only if c > c* = 2a/ /'(0). If we look at the equation 

u t = DAu+f{u), (1.4) 

for some positive constant D, then a simple rescaling yields that the minimal KPP speed 
of (ll.4p is given by c* D = 2\JT)^J f'(0). One can then notice that, in a homogeneous setting, 
the minimal speed is increasing with respect to the diffusion factor (the map D i— > c* D is 
increasing over (0, +oo)). 

1.0.2 Heterogeneous framework, notations and setting 

The result of Kolmogorov et al. [6] has been generalized to media with periodic spatially 
dependent coefficients (see p], [H], [13], [14] for example) and to settings with space-time 
dependent coefficients (see [2], [9j [10] ) . We recall here some results which are very relevant to 
this present work and we also introduce the setting that will be considered. The equation 
that we will consider is of the type 

ut(t,z) = V z -(A(z)V z u) + q(z)-V g u + f{u), teR,zett, (1.5) 
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where the domain Q : = R x u is an unbounded C 3 open connected subset of M. N =lxK 
(N > 1). We represent each z G Q as 



In order to cover all possible cases, we will now give a generic description of uj which 
appears in Q when the space dimension is iV > 1. The set u is assumed to have the 
following structure. There exists an integer d G {0, . . . , N — 1} and L\ > 0, . . . , Ld > so 
that an element y G u is represented by (2/1, 2/2) G R d x R^ -1 ^ and 



where (ej)i<j<Ar is the canonical basis of M N . 

• Notice that the case d = corresponds to u cj is bounded" and hence Q = R x u; is 
an infinite cylinder of section u. 

• In the cases where d > 1, the set w is unbounded. 

As particular examples of f2, one can have the whole space R N = R x R 7 ^ 1 in which case 
d = N — 1. One can also have the whole space R^ except a periodic array of holes (periodic 
perforations). In such case, d is also equal to AT — 1. For N = 2, one has d G {0, 1}. The 
case where d = 1 (with N = 2) means that u is unbounded connected and satisfies (II. 6p 
and thus Q = R 2 . For N = 3, d can be 0, 1, or 2. The case where d = 1 corresponds 
to u C R 2 bounded in one direction and unbounded in the other one. In the unbounded 
direction, u has to be periodic with a period that we denoted by Li in the general setting 
above. The case where d = 2 corresponds to d — N — 1, and thus, has been discussed 
above. 

We mention that in all cases (even when d — 0), the domain f2 := R x 00 has a 
periodicity cell which we denote by C := [0, 1] x C w where 



{C w = oj, when d = 0, 
C w = {(yi, 3/2) = (Vi, • • • , 2/i, 2/2) G w, 2/1 G [0, Li] . . . ,yf E [0, L d ]} otherwise. 

Concerning the reaction term, we will only deal with nonlinearities / depending on 
u in order achieve the proof of non-monotonicity of the speed with respect to diffusion. 
Precisely, / is of the homogeneous KPP type f ll.2tll.3p . 

The advection term, in this paper, is assumed to be a shear flow. That is a vector field 
q(x, y) = (qi(y), 0, . . . , 0) of class C 1,a (f2) (a > 0). The advection is assumed to satisfy 



z = (x, y), x G R, !/6w. 




(1.6) 




. . . , Ld)— periodic with respect to y (whenever d > 1) 



(1.8) 
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Obviously the above vector field satisfies 



V^.y • q = in Q and q ■ v = on dQ, 

where u is the outward normal on dfl and it is given by v = (0, z/ w ) where u u is the outward 
normal on du. 

Let us now describe a wide class of diffusion coefficients for which the existence of 
traveling fronts and minimal speed of propagation holds according to [T]. The diffusion 
matrix A(x,y) = A(x,yi,y 2 ) = (Aij(x,y))i<ij<N is a symmetric C 2 ' a (fl) (with a > 0) 
matrix field satisfying 

A is (1, Li, ■ ■ ■ , Ld)— periodic with respect to (x, y±), 
30<a 1 <« 2 ;V(x,|/)Gfi 1 VeGl Ar , (1 g) 

l<i,j<N 

In the above setting where Q = R x u satisfies (jl.6p . q is of the type (11.81) and A 
satisfies (11.91) . one can rewrite (11.51) as 

u t (t,x,y) = V • (A(x,y)Vu) + q x {y)d x u + f(u) inlxKxu, (1.10) 

together with the boundary condition (in the cases where duo ^ 0) 

u(x, y) ■ AV XjV u(t, x, y) = 0, (t, x, y) <E R x R x du. (1.11) 

Now in this non-homogenous setting, we set e = (1, 0, . . . , 0) G M N as the direction of 
propagation. We recall, from [TJ [131 EL that a pulsating traveling front in the direction of 
— e that propagates with a speed c G R is a classical solution u = u(t, x, y) := (p(x + ct, x, y) 
that connects to 1 as follows 

lim u(t, x, y) = and lim u(t,x,y) — l 

x—¥—oo x— >+oo 

(locally in t and uniformly in (x, y)), and satisfies u{t H — , x, y\ + k, y 2 ) = u(t, x + 1, yi, y 2 ) 

for any k G L{L - ■ ■ x L^Z. In terms of 0, the latter means that is 1-periodic in x and 
(Li, • • • , Lrf)— periodic in y\. Under the assumptions (ll.6H1.81fL~9l) on Q, q and A, and having 
a KPP nonlinearity ( ll.2til.3p one knows from [lj that there exists a minimal speed c* : = 
c *A,n,qj( e ) > so that a pulsating traveling front with a speed c exists if and only if 
c > c* n A q j{e). The minimal speed c* has been expressed in terms of the coefficients of the 
reaction-advection-diflusion problem via a variational formula which involves eigenvalue 
problems in j3j and [13] . This formula has been used in many works to study various 
asymptotic and homogenization regimes of pulsating traveling fronts (see for eg. [1], [5], 
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2 The non-monotonicity of the minimal speed with 
respect to diffusion (in the presence of an advection 
field) 

Berestycki, Hamel, and Nadirashvili [3] proved (in part 2 of Theorem 1.10) that, having 
any periodic domain Q C M. N , q = and a constant growth rate f' u {x, y, 0) (which holds in 
the case / = /(it)), the map (3 i— > Cq p AO j(e) is increasing in (3 > 0. Nadin [8] has proved, 
in Theorem 2.5, that the same result still holds even when f' u {x, y, 0) depends on the spatial 
variables (x,y). However, we notice in both results of J3J |8] the absence of an advection 
term. 

Furthermore, in Theorem 6.1 of [4J, where Q = M x u, q is a shear flow of the form 
q(x, y) = (qi(y), 0, . . . , 0) in Q, while A and / satisfy certain assumptions which are valid 

in our present setting, it was proved that the map j3 i y ^= is decreasing with 

V P 

respect to /3 > in both cases qi ^ or qi = over u. Having those results in [3], [1] and 
[8], there arise naturally the following two questions. 

• First: Do we still have the increasing behavior of the minimal speed with respect 
to the diffusion factor /3 in the presence of an advection, even if the nonlinearity is 
homogenous? 

• Second: Owingto Theorem 1.1 in [3J, the map A i-> /( e )' where A varies in the 
ordered family \j of positive definite matrices satisfying ( |L9p is well defined (provided 
that Q, q and / satisfy (11. 6p . (II. 8p and (ll.2tll.3p ). We investigate the variation of the 
minimal speed of propagation with respect to the matrix of diffusion. More precisely, 
if A and B are two positive definite matrices satisfying (ll.9p and if A < B, do we 
still have £& iVW (e) < c* UBtq j{e) whenever q £ 0? 

As a matter of fact, the presence of a shear flow will change the monotone behavior of 
the minimal speed with respect to diffusion. We prove in Subsection 12.11 that the answer 
is negative in general even when the diffusion matrices A and B are proportional to the 
identity matrix and the nonlinearity is homogenous. We give a counterexample when the 
advection is large (up to a scaling, this gives a counterexample to the monotonicity of the 
speed with respect to diffusion with a prefixed term and a small reaction). In Subsection 
12. 2[ we prove that, for a fixed nonzero shear flow and a fixed reaction, the answer to the 
second question is negative in general for matrices A < B which are not equal up to a 
positive scalar. We mention that the second result cannot be seen as a trivial consequence 
of the first one by using the argument that "a small perturbation of the strict inequality 
(12. 2 p remains a strict inequality and one can then perturb the diffusion matrices" (see 
the more precise details at the beginning of Subsection 12. 2p . Lastly, we mention that our 
counterexamples apply when the diffusion matrices do not depend on the variable x. One 

1 We say that A — A(x,y) < B — B(x,y) if and only if for each (x,y) £ f2 and for each £ £ M. N , we 
have £ • A(x, y)£ < £ • B{x,y)S > . Also, we say that A < B if and only if for each (x,y) £ and for each 
£ £ R N \ {0}, we have £ • A(x, y)£ < £ • B(x, y)£. 
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knows that, in such case, the traveling fronts of the reaction-advection-diffusion problem 
fll.lOP will have the form <p(x + ct, y) and they will be invariant with respect to the frame 
moving in the direction — x (or — e). Moreover, we know that in a homogenous setting the 
speed is monotone with respect to diffusion factors (see Section [1.0. ip . Thus, the closer we 
are to a homogeneous setting the harder it is to prove non-monotonicity of the KPP speed 
with respect to diffusion. Indeed, the setting where we construct our counterexamples can 
be taken very close to the homogeneous one with only one heterogeneity coming from the 
advection. This gives sharpness to our results. 



2.1 A counterexample devoted to answer the first question 

In this subsection, we will show a reaction-advection-diffusion problem whose diffusion 
matrix varies in the subfamily of positive definite matrices PDj^ = {(3 Id, (3 > 0} (where 
Id stands for the N x N identity matrix) while a shear flow will make the Part 2 of 
Theorem 1.10 [3j not valid anymore, even if the nonlinearity / is taken homogeneous. 
Let's first announce the following result. 

Theorem 2.1 Let e = (1, 0, . . . , 0) G R N , tt = R X u C R N satisfying / TOP where u may 
or may not be bounded. Assume that the nonlinearity f = f(u) is a homogenous "KPP" 
nonlinearity satisfying U.£M1.3\) and let's start with a shear flow q(x,y) = (qi(y),0, ■ ■ ■ , 0) 
defined over Q such that q\ ^ (satisfying U.8\) when d > 1) and so that 

< 2y/f'(0) + 5 < max (qi{y)) - 5, for some S > 0. (2.1) 

Then, there exists M x := M 1 (5) > for which there corresponds < E\ := Sx(8, Mi) < Mi 
such that 

V0<e< £l <M 1 , VMv^ 9 ,/ (e)> Wi^,/ (e) - (2 ' 2) 

Proof. Before going further into details, we recall that the variational formula of the KPP 
minimal speed (see [3] for example) yields the continuity of (k, p, fi) kA af( e ) w hh 

respect k > 0, \i > and pel. This continuity with respect to the factors of the reaction, 
diffusion and advection will be useful to construct our proof. Due to the presence of several 
parameters in the problem, we proceed in simple steps. 

Step 1. In the author's work [4], an asymptotic regime for the speed c* within large 
diffusions MA together with advection fields of the form M 7 g was proved (where M — y 
+oo). The exponent 7 was allowed to be any number in the interval (0, 1/2]. Precisely, 
Theorem 4.1 of [1], with 7 = 1/2, yields that 

c* ^- (e) 

fl,MId,y/Mq, f 



lim "'^"'^^ =2 vTm 

Using the above together with the continuity of c* with respect to M, there then exists 
M := M {8) > such that 
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Step 2. We fix M 1 > max(l, M (S)). Then, 

< c k Mi u, jw iq , ,<.<)< 



(2.3) 



Step 3. For the fixed number Mi, we also have a/Mi q in the form of shear flow. We 
now look at the effect of small diffusion while y/M[ q is considered as an advection field. 
Indeed, Theorem 3.3 of El Smaily jl] yields that 

lin ! c O p ia rw „ f ( e ) = ^^(VMiqi(y)) = \[M~\ max(gi(>)). 
£ -s-o+ ll,£ la, y M\ q, j yeui y eu 

Consequently, there exists E\ = £i(S, Mi) G (0, Mi) (we can choose e\ as small as we like) 
such that 



> a/mT 

Step 4. Owing to (J2HJ, Q and (12^1) . one then gets 



max(gi(y)) - 5 



> 0. 



(2.4) 



V0<£<£ b c* , rj-j- Je)>c* , rj-r .(e) 

11, e Id, a/Mx g, / v y \l, Mild, \J Miq, f K ' 



and this completes the proof. 



□ 



2.2 Case of non-proportional diffusions (answer to the second 
question) 



In Subsection 12.11 we saw that a prefixed advection field was multiplied by a large enough 
constant Mi in order to spoil the increasing behavior of the minimal speed with respect to 
diffusion matrices which are proportional to the identity matrix. A simple scaling of the 
reaction-advection-diffusion equations that correspond to Theorem 12. II implies that, under 
the same notations, 

°h, (e/VM^Id, q, f/y/Ml (e) > °h, JWjd, q, f/sfM~S e) ' 

and by continuity of the minimal speed with respect to diffusion coefficients, one can then 
find non-proportional matrices A and B (in the neighborhood of (s/ \fM[)Id and yM\Id 
respectively) that satisfy the general assumptions of Section [TJ and such that A < B and 

C h, A, q, f/^Ml {e) > °h, B, q, f/y/M^' 

This can be summarized as a perturbation argument applied to a strict inequality. However, 
we notice that the latter inequality corresponds to reaction terms which are small. The 
importance of the result we give in this subsection is that it leads to the construction of 
a counterexample to the monotonicity of the speed with respect to diffusion while q and 
/ are prefixed (neither small nor large). It turns out that the counterexample given in 
Theorem 12.31 concerns diffusion matrices which are non-proportional by construction and 
not by a perturbation argument of the result of the previous section. 
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Notation 2.2 For each b > 0, let A b denote the N x N symmetric positive definite matrix 
having the form 





(10 


o \ 




b '■■ 




A b = 















\ 


b J 


Theorem 2.3 Let e = (1, 0, . . . , 0) 


GR N ,0 = fxw 


C R N 


or unbounded satisfying $1.6\), and 


let q = (qi(y),0, ■ 


..,o) 



where u can be either bounded 



where qi ^ on u. Assume that the nonlinearity satisfies U.2W1.3\) . For each b > 0, 
consider the reaction- advection- diffusion problem 



Ut{t,x,y) = V ■ (A b Vu) + q 1 {y)d x u + /(«), 

= d xx u + bAyU + qi(y)d x u + f{u) inM. x Q, 
u(x, y) ■ A b V x , y u = Vu{y) ■ V y u(t, x, y) = for (t, x, y) G R x R x du, 

where A b is the matrix introduced in Notation \2.2\ Then, 

1. 



lim c* a A 

b^+oo ^ l ^b,qiJ 



2. 



lim c* 



'T + u n A h a f ( e ) = m -^ x (?i + 2 v/fM- 
Proof. Consider the following change of variables 



x 



V(t, x, y) G R x R x a, v(t, x, y) = u(t, y). 
Thus, V(t, x, y) G R x R x w, 

v t (t, x, y) = u t (t, -^=, y), d x v(t, x, y) = -j=d x u(t, -^=, y), 

1 X X 

d xx v(t, x, y) = -d xx u(t, -j=, y) and A y v(t, x, y) = A y u(t, — , y). 

Owing to the invariance of Q in the x— direction, we have the boundary condition 

V(t, x, y) G R x dtt, u(x, y) ■ V XiV v(t, x, y) = 0. 

Consequently, the problem ( 12.51) is equivalent to the problem 

v t (t,x,y) = bd xx v + bA y v + \fbq x (y)d x v(t,x,y) + f(v), 
= bA XtV v + Vbqi(y)d x v + /(f) in R x R x u, 
v(x, y) ■ V XjV v(t, x,y) =0 for (t,x,y) G R x R x du. 



(2.5) 



(2.6) 
(2.7) 



(2.8) 



The minimal speed of problem (12. 8p exists in this setting (according to Theorem 1.2 

in [1]) and is denoted here by c* /- (e), where 7d stands for the identity N x N 

it, bid, y og, j 

matrix. As in the proof of the Theorem 12.11 we apply Theorem 4.1 of jl] with 7 = 1/2 to 
conclude that 

c* (e) 
lim Wl^kil = 2Jm. (2.9) 

On the other hand, looking at the problems ( 12. 5ft and ( 12. 8ft and the relation between w 
and f and owing to the minimality (uniqueness) of the threshold c*, one knows that the 
corresponding minimal speeds have the following relation 

Vt>0 ' Wv^,/ (e) = V H^ 9 ,/M- (2-10) 
Together with (12.91) . we obtain that 



Proof of ( 12. TP . For the limit as 6 — > + , we know from the rescaling formula (I2.10p 
that proving (12. 7| ) is equivalent to prove that 

c* (e) 
bid, yfbq, f 



ii, ui u, yu q, = ^ + 2 ^^y^ 

6-»0+ a/5 w 

Indeed one can apply similar techniques to those used in the proof of Theorem 3.3 of jl] 

(C replaced by /'(0)) to obtain prove this. However, for the sake of completeness, we do 

the proof here in details. One knows that c* r (e) is given by the variational 

\l, bid, Vbq, f 

formula of [3] 

V6>0 ' c n uta ru ,(e) = min^p, (2.12) 
fi, 6/<i, Vbq, f \>o X 

where k(X) is the principal eigenvalue of the following problem 

Ln, e ,bid, q j,x 4> = bAtp + 2b\e ■ + Vb~q 1 {y)d x ip 

for- 1 (2-13) 



+ 



^ in f2 = R x 



with the boundary conditions v ■ Vtp = \v ■ eip on dQ. The eigenfunction ip = ip(x,yi,y2) 
is positive in fl, 1-periodic in x, and (L 1; . . . , L d )— periodic in y 1; and unique up to mul- 
tiplication by a non-zero constant. Now, having the coefficient qi independent of x, this 
yields that one can choose ip independent of x. Hence, the elliptic operator L n e bId ^ q . x 
can be reduced to the symmetric operator 
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Consequently, we have the following Rayleigh quotient over if^C^), where C u is the peri- 
odicity cell of u introduced in Section [1] above, 



VA > 0, V6 > 0, 



-k(X) = min 



bf \V V \ 2 dy - XV% I qM^-l [\ 2 b + f'(0)]v 2 (y)dy 



<p {y)dy 



Formula (I2.14p leads directly to the upper bound 



(2.14) 



VA>0,V6>0, M^< ma _ x (g 1 ( 2/ )) + Av / 6 + {^ 
XV b XV b 



Testing the right hand side at A = \ ^-4^ , one then concludes that 



.:. MA) 



V6 > 0, mm < ma_x(g 1 ( 2/ )) + 2y/f(fi), 
A >° XVb 

and hence 

c* (e) 
limsup to,Md,y/bq, f — < max(g 1 (y)) + 2^/fiO). (2.15) 

6^o+ vb y euj 

For the sharp lower bound of the liminf as b — » + , we proceed as follows. Let yo be a 
point of C cJ where the periodic function gi(y) attains its positive maximum (we recall 
here that q\ ^ 0, q± is periodic in y and has a zero average over C w and thus the point yo is 
in the interior of the cell CJ). Then, by the definition of a maximum point of a continuous 
function, for each < 5 < qi(yo), there exists a neighborhood Us C C w of y sucn that 

VyeU s , qi(y)>qi(y o )-S>0. 

Now we test the Rayleigh quotient against a smooth function <f)(y) compactly supported 
in Us and normalized by j c <fi 2 dy = 1 to get (after switching from — A;(A) to +k(X)) 

^ > Xb + 1/3(6) + Vb( gi (y ) - 6), (2.16) 

where 0(b) = /'(0) - 6 / | V0| 2 . We notice that (3(b) > for < b < b := / ^ ■ The 

Jc/ * / |V0| 2 

function A i — ^ A6 + T/^(&) attains its minimum over (0, +oo) at A* = yj and the value 
of this minimum is 

2V^vW)- Referring to fl2~T6l) . one then gets 

V0<6<6 ,VA>0, ^>2Vb^/W) + ^b(qi(yo)-s). 
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The right hand side of the last inequality is independent of A. Thus, the variational formula 
( jgZE2D yields that 

c* (e) 
bid, y/bq, f 



V0 <6<6 , »^ >2VW + (gi(y )-g), 

for all < 5 < qx{yo). Thanks to the non-dependence of Ug on 6, we pass to the limit as 
b ->■ 0+ to get 

c* r (e) 

V0 < 5 < ftdft,), liminf ^MMMil > 2v /f(0) + (gi(y ) - 5). 

This is enough to conclude that 

c* (e) 
fi, bid, Vbq, f 



hminf "' Ma ' V - 0g ' J > 2V? 7 (0) + (ft(ito)) (2.17) 
and, together with f !2. 15j) . finish the proof of (12. 7p . □ 

A counterexample as an application of Theorem 12.31 

Let e = (l,0,...,0) and H = 1 x w (in particular Q can be taken as the whole space M ). 

Choose / = f(u), and q = (q 1 (y),0, ■ ■ ■ ,0) with / q 1 (y)dy = and q 1 ^ 0. Thus there 

Jc 

exists S > such that 



2 % /7(0) + ^<niax(g 1 (y)) + 2V7 7 (0)-5. 

yew 

It follows, from Theorem 12.31 above, and from the continuity of the speed with respect to 
the diffusion factor (variational formula of [3]) that there exist 6q > and M > such 
that 

\f0<e<e , c* A f (e) >ma_x(g 1 (y)) + 2 A /7 7 (0)-5 and 

VM>M >0, c* nAMjqJ (e)<2^rW) + S. 

Consequently, choosing e small enough and M large enough, it follows that A M > A e 
in the sense of order relation on positive definite matrices; however, 

C kA M ,qJ^< C kA £ ,qJ^- 

Therefore the answer to the second question is negative in general even when the 
advection is a fixed shear flow and the nonlinearity / is fixed and homogenous. □ 
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